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The doublet channel neutron-deuteron scattering amplitude is calculated 
in leading order effective field theory (EFT). It is shown that this amplitude 

t^J- \ does not depend on a constant contact interaction three-body force. Satis- 

factory agreement with available data is obtained when only two-body forces 

q ' are included. 
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An intriguing difficulty arises in the application of leading order EFT to the three-body 



problem. One finds that the full amplitude describing three-boson scattering, or nucleon- 
deuteron (nd) scattering in the J = 1/2 channel, is sensitive to the cutoff used to solve the 
scattering equations - even though each perturbation diagram, with resummed two-body 
interactions, is individually finite. In |l[] it was argued that the addition of a one-parameter 
three-body force counter-term is necessary and sufficient to eliminate this cutoff dependence. 
On the other hand, in refs. || and @] we have shown, on the example of three bosons, that the 
cutoff dependence is just a natural consequence of the existence of infinitely many solutions 
to the given scattering equation; moreover, by carefully identifying the physical amplitude 
from amongst the infinitely many non-physical solutions, we have shown that the cutoff 
problem can be solved without the introduction of a three-body force. 

In the present contribution we show that the physical three-body scattering amplitude 
in fact does not depend on the constant contact interaction three-body force at all. Further, 
we demonstrate that for doublet channel nd scattering, good agreement with experiment is 
obtained with the inclusion of two-body forces only. 



II. WHY THERE IS NO THREE-BODY FORCE DEPENDENCE 

To show that the leading order EFT three-body amplitude is independent of constant 
three-body forces, it is sufficient to restrict the discussion to the case of three bosons. 

A. Scattering amplitude without three-body forces 

In the three-boson case without three-body forces, the s-wave particle-bound-state scat- 
tering amplitude a(p, k) satisfies the equation 

a(p,k) = M(p,k) + — dqM(p,q)^-f- -o(g,fc), (1) 
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where 



M(p,g) = ~ ll + ^p2-m£ 



1 i ( q + pq + p — mE 

2pq \q 2 — qp + p 2 — mE 



(2) 



In this equation k (p) is the incoming (outgoing) momentum magnitude, E = 3k 2 /4m — 
1 /ma 2 , is the total energy, and 02 is the two-body scattering length. Here it is assumed that 
the summation of perturbation theory diagrams and loop integration can be interchanged 
in the sense that the difference is of higher order and hence negligable in given leading order 
calculations. In general such assumptions have to be investigated very carefully as they may 
lead to fictitious fundamental problems |2|. 

Eq. ([[]) is known as the S-TM equation ||, and in the three-boson case has A = 1. Three 
nucleons in the spin J = 1/2 channel obey a pair of integral equations with similar properties 
to this bosonic equation, while the J = 3/2 channel corresponds to A = —1/2. For A > 
Danilov's work [7[] shows that the homogeneous equation corresponding to Eq. (|I|) has a 
solution for arbitrary E; in particular, there exists a solution for every energy corresponding 
to the scattering of a projectile off a two-body bound state. 

The existence of these solutions implies that Eq. ([|) has an infinite number of solutions. 
In fact the homogeneous equation has more than one solution for any given E. Writing these 
solutions as a\ where % — 1, 2, 3, . . ., the most general solution of Eq. ([!]) can be written as 
a = a p + J2i Cia\ where a p is any particular solution. It is useful to examine the asymptotic 
behaviour of a(p, k) for large p. Because the inhomogeneous term M behaves asymptotically 
as 1/p, it follows that either (i) a — > faster than 1/p, or (ii) the asymptotic behaviour of 
a is determined by the asymptotic behaviour of the homogeneous solution a^. In the latter 
case the asymptotic behaviour has the form 

a(p,k)=J2A(k)p s * + 0(l/p) (3) 
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where Si are roots of the equation 

x _8A sinW6 =Q (4) 
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The summation in Eq. (|3|) goes over all solutions of Eq. (Q) for which |Res| < 1. For A = 1 
Eq. (|^) has two roots for which |Res| < 1: s — ±is , where s ~ 1.00624, so that Eq. @ 
gives the asymptotic behaviour of the amplitude as 

a(p,k)^A 1 (k)p iS0 + A 2 (k)p- iS0 . (5) 

By contrast, for A < the homogeneous equation has no non-trivial solution and the solution 
of Eq. ([[]) is unique; in this case the physical amplitude a must vanish asymptotically faster 
than I /p. 

In refs. |3],|J we have shown that the oscillatory behaviour of the general amplitude a for 
A = 1 is simply an artifact of the homogeneous equation (corresponding to Eq. (p])) having 
non-zero solutions, and that the physical amplitude does not display this spurious behaviour; 
instead, it behaves just like the solution for A < 0, namely, it vanishes asymptotically faster 
than 1/p. Thus amongst the solutions given by Eq. ([5]), the physical solution is the one with 
A\(k) = A 2 (k) = 0. More generally, for any A > the physical amplitude is the one that 
has no admixtures of homogeneous equation solutions, and by the above argument, it must 
therefore vanish asymptotically faster than 1/p. 

Considering Eq. (HD for the case where a is the physical amplitude, since the free term of 
Eq. ([I]) behaves like 1/p for large p, the coefficient of 1/p coming from this inhomogeneous 
term should cancel the coefficient of a similar term coming from the integral part (we note 
that this argument is valid for both A = 1 and A = —1/2). Hence 

= —^ + —^ -i, pj a(q,k). (6) 
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B. Identical scattering amplitude with a constant three-body force 



Multiplying Eq. (|) by 2/y/3(l/a 2 + fi/Ap* - mE)H and adding the result to Eq. (|), 
we obtain the scattering equation where the constant three-body force H is included to all 
orders (an H is simply added to the term in the square bracket of Eq. (^j). Hence the physical 
non-oscillating solution of Eq. ([!]) with no three-body force also satisfies the modified Eq. ([!]) 
where an arbitrary H is included. Hence the inclusion of a constant three-body force has 
no effect on the physical scattering amplitude. 

In a recent paper || it has been shown that doublet channel neutron-deutron scattering 
amplitude in EFT with effective range parameters taken into account, does not exhibit any 
dependence on a constant three-body force. This result is therefore in agreement with the 
observations of the present work. 

III. ND SCATTERING WITH TWO-BODY FORCES ONLY 

Doublet channel neutron-deuteron scattering in leading order EFT is analogous to the 
scalar case and does not involve any additional problems. Starting from the EFT Lagrangian, 
one can obtain the following system of doublet channel neutron-deuteron scattering equa- 
tions 
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FIG. 1. Doublet channel nd scattering phase shifts. Data are from phase shift analysis of van 
Oers and Seagrave (dots) and a measurement by Dilg et al. (triangle). 



2 r°° 
a(p, k) = M(p, k) + - dq M(p, q) [G 1 (q)a(q, k) + 3G 2 (q)b(q, k)] , 
TV JO 

2 r°° 
b(p, k) = 3 M(p, k) + - dq M(p, q) [3Gd(g)a(g, k) + G 2 (q)b(q, k)} , 

7T J0 

where a and b are the neutron- 3 ^ (nd) and nucleon- 1 5o amplitudes, respectively, 
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are the two-particle scattering lengths in singlet and triplet channels, and E = 3k 2 /4m — 
l/m{a t 2 ) 2 is the total energy. Apart from a factor of 1/4, the driving term M(p, q) differs 
from Eq. (0) only in that a 2 is replaced by a\. 

Solving the system ©-© and isolating the non-oscillating solution for a(p, k) we obtain 
the physical amplitude for nd scattering. The results are shown in Fig. 1. In order to isolate 
the physical amplitude we had to solve the corresponding homogeneous system of equations. 
This task is especially difficult due to the existence of a continuum of solutions corresponding 
to the continuous spectrum of scattering energies. The method employed to achieve this 
numerical solution does not allow us to obtain high accuracy for low momenta ||^]] - that 
is why our curve does not extend to the origin. We note that while our calculations fit the 



experimental data quite well, the accuracy of these data is open to question. Ref. |10[ does 
not contain error estimates and ref. [TTJ claims that at least the scattering length calculated 



in ref. [TTJ] may be incorrect. 
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